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This  Technical  Memorandum  presents  the  results  of  a maximum  likelihood 
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catapult  launch  configuration.  The  aerodynamic  coefficients  presented  herein  were 
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the  EA-6B  airplane. 
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INTRODUCTION 

BACKGROUND 

1.  NAVAIRTESTCEN  was  directed  by  reference  I to  conduct  flight  tests  to 
determine  the  single  engine  catapul*  launch  capabilities  of  the  EA-6B  airplane. 
Tests  were  conducted  to  determine  the  single  engine  gross  weight/temperature 
catapult  laimch  envelope,  the  optimum  piloting  techniques  and  procedures  to  be 
utilized  in  the  event  of  cm  engine  failure  during  or  immediately  following  catapult 
launch,  and  the  feasibility  of  using  20  rather  than  30  deg  flaps  for  launches  at 
heavy  weights  and  high  ambient  temperatures.  The  results  of  these  tests  were 
reported  in  reference  2.  In  addition,  tests  were  authorized  to  determine  the  EA-6B 
airplane's  pitching  moment  coefficient  in  order  to  further  study  the  single  engine 
catapult  launch  environment  through  the  use  of  computer  simulation.  This 
Technical  Memorandum  presents  the  normal  and  2ucial  force,  as  well  as  the  pitching 
moment  and  stability  and  control  derivatives  obtained  from  a maximum  likelihood 
analysis  of  flight  test  data  in  the  catapult  launch  configuration. 

FUGHT  TEST  DATA 

2.  The  test  data  were  obtained  during  two  separate  flights  at  an  average  altitude 
of  7,000  ft  (2  130  m).  One  flight  was  in  the  clean  (no  stores)  loading  eind  one  with 
three  tactical  jamming  (TJ)  pods  and  two  300-gal  (1  135.6 S,)  external  fuel  tjmks 
installed.  The  test  airplane  utilized  in  this  program,  BuNo  156478,  was 
aerodynamically  representative  of  a production  EA-6B  airplane  with  the  exception 
of  the  installation  of  a test  nose  boom  containing  an  NACA  pitot-static  head,  Kiel 
total  pressure  probe,  and  angle  of  attack  and  angle  of  sideslip  vanes.  The  geometric 
and  moment  of  inertia  characteristics  of  the  test  airplame  in  the  catapult  launch 
configuration  are  presented  in  appendix  A. 

3.  The  evaluation  was  conducted  in  an  aerodynamic  configuration  representative 
of  the  catapult  laimch  environment  (landing  gear  DOWN,  flaps  30  deg,  speedbrake 
IN).  Thrust  for  level  flight  was  utilized  rather  than  the  catapult  launch  setting  of 
Military  power  in  order  to  facilitate  establishing  trimmed  level  flight  conditions. 
This  did  not  have  a significant  effect  on  the  estimated  parameters.  However,  the 
test  stabilizer  trim  requirements  were  slightly  less  than  those  which  would  be 
obtained  at  full  Military  power  due  to  the  reduced  downwash  characteristics  at  the 
horizontal  stabilizer.  The  stabilizer  required  to  trim  with  thrust  for  level  flight  is 
presented  in  figure  1. 

h 

i: 


1)  NAVAIR  msg  152018Z  May  1975 

2)  NAVAIRTESTCEN  Technical  Report  SA-105R-76,  Evaluation  of  Single  Engine 
Catapult  Launch  Capabilities  of  the  EA-6B  Airplane,  Final  Report,  of 
10  February  1977 
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ALTITUDE -7J>00  FT  (2  130  m) 

CG  RANGE  - 23^  TO  Zi2  % MAC 
LANDING  GEAR  - DOWN 
FLAPS -30  DEG 
THRUST  FOR  LEVEL  FLIGHT 
• O 3TJ  PODS  AND  2 DROP  TANKS 
■ □ CLEAN  LOADING 


Lb 


SASOFF 

SASON 


PDWER  OFF  LIFT  COEFFICIENT 


Figure  1 

Level  Flight  Trim  Requirements 

4.  The  test  airplane  was  instrumented  for  longitudinal  flying  qualities  data  as 
listed  in  appendix  C.  Test  data  were  telemetered  to  the  NAVAIRTESTCEN  Real- 
Time  Telemetry  Processing  System  and  recorded  on  onboard  mangetic  tape.  The 
pilot  techniques  required  to  generate  the  flight  test  data  necessary  for  the 
parameter  identification  process  were  standard  dynamic  longitudin2d  stability  test 
techniques.  The  stability  2md  control  derivatives  describing  the  short  period  mode 
of  motion  were  obtained  from  the  response  to  pitch  doublets  and  sinusoidal  stick 
pumping  inputs.  The  stability  derivatives  describing  the  variation  of  forces  and 
moments  with  respect  to  forward  velocity  were  obtained  from  stamdard  phugoid 
maneuvers.  All  tests  were  conducted  from  stabilized  trim  flight  conditions. 
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PARAMETER  IDENTIFICATION  ALGORITHM 

5.  The  parameter  identification  technique  used  in  this  analysis  is  a maximum 
likelihood  identification  algorithm  (SCIDNT)  developed  in  a joint  program  by 
NAVAIR,  Office  of  Naval  Rese2U'ch,  NAVAIRTESTCEN,  and  Systems  Control, 
Incorporated  (SCI).  Maximum  likelihood  estimation  theory  is  discussed  in  references 
3 and  4 and  details  of  the  computer  cdgorithm  used  in  this  analysis  are  contained  in 
reference  5.  This  maximum  likelihood  algorithm,  ais  implemented,  is  basically  a 
linear  program  capable  of  estimating  linear  airframe  stability  and  control 
derivatives.  The  maximum  likelihood  identification  procedures  used  provide  for 
simultaneous  estimation  of  aerodynamic  coefficients  and  instrumentation  errors. 
Thus,  predata  filtering  or  correction  for  instrumentation  errors  is  not  required  prior 
to  the  parameter  estimation  procedure.  The  short  period  and  phugoid  response 
characteristics  were  analyzed  at  25  and  5 samples  per  second,  respectively. 

6.  The  equations  of  motion  (body  fixed  axis  system)  used  in  this  analysis  are 
presented  in  appendix  B.  The  body  oucis  system  used  is  defined  as  follows: 

a.  Orthogonal,  right-hand  system,  fixed  to  aircraft. 

b.  Origin  fixed  at  aircraft  eg. 

c.  X-axis  is  in  the  plane  of  symmetry  fixed  to  fuselage  reference  line,  positive 
forwcurd. 

d.  Y-axis  is  perpendicular  to  the  plane  of  symmetry,  positive  to  the  right 
(right  wing). 

e.  Z-axis  is  in  the  plane  of  symmetry,  perpendiculcu:  to  the  X-axis,  positive 
downwcurd. 

7.  The  final  set  of  aerodynamic  coefficients  presented  in  this  report  were 
determined  by  combining  the  results  of  numerous  memeuvers  at  the  same  flight 
condition.  The  parameter  estimates  will  vary  between  maneuvers  because  of  the 
statistical  nature  of  the  random  errors  that  occur  in  the  measurement  system  and 
the  different  amounts  of  information  in  each  mcuieuver  concerning  the  peurameters 
to  be  estimated  (parameter  identifiability).  This  vau-iation  in  information  or 
parameter  identifiability  is  related  to  the  degree  of  excitation  of  the  modes  of  a 
system  by  a particulcu:  input  and  the  corresponding  ability  to  accurately  identify 
the  parameters  of  the  system.  In  other  words,  in  order  to  identify  a mode  (or 
parameters)  of  a system,  that  particular  mode  must  be  excited  by  the  input. 


3)  NAVAIRTESTCEN  Technical  Report  FT-77R-73,  First  Interim  Report,  Devel- 
opment of  Digital  Airframe  Parameter  Identification  Technology,  of  21 
January  1974 

4)  NASA  Contractor  Report  NASA  CR-2200,  Maximum  Likelihood  Identification 
and  Optimal  Input  Design  for  Identifying  Aircraft  Stability  and  Control 
Derivatives,  of  March  1973 

5)  Systems  Control,  Incorporated,  Engineering  Report,  SCIDNT  I Theory  and 
Applications  - Technical  Report  No.  3,  of  December  1974 
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Parameter  estimates  will  vary  from  rim  to  run  depending  on  the  repeatability  of 
the  input  or  types  of  inputs  used.  Thus,  eis  outlined  in  reference  6,  the  estimates 
from  the  different  maneuvers  must  be  combined  to  make  full  use  of  available  data. 

8.  In  terms  of  statistical  characteristics,  parameter  estimate  accuracy  is 
discussed  in  terms  of  confidence  bounds  and  estimated  parameter  error  (EPE). 
Parameter  estimate  confidence  bounds  cure  determined  in  the  identification 
algorithm  based  on  the  Cramer-Rao  lower  bound  (references  4 amd  5).  A la 
confidence  bound  represents  a 67%  confidence  interval  statement  of  the  parameter 
estimates  accuracy,  and  a 2a  confidence  bound  is  a 95%  confidence  interval.  For 
example,  a 2a  confidence  bound  implies  that,  in  a large  number  of  flight  tests 
under  identical  conditions,  95%  of  the  time  the  parameter  estimate  will  be  within 
+ 2a  of  the  true  parameter  value.  Parameter  estimate  confidence  bounds  are  used 
as  a guide  in  assessing  the  accuracy  of  individual  parameter  estimates;  however, 
this  confidence  bound  does  not  predict  anything  about  the  outcome  of  a single 
experiment  in  a definite  sense.  Thus,  in  addition  to  parameter  estimate  confidence 
bounds,  other  criteria  must  be  considered,  such  as  static  stability,  time  history 
comparisons  (response  prediction  criteria),  amd  modal  response  frequency  and 
damping  calculations  in  assessing  parameter  estimate  accuracy.  The  relative 
accuracy  of  all  parameter  estimates  is  qualitatively  assessed  in  this  report  by  using 
the  EPE  value  which  is  defined  as: 


EPE  (Estimated  Parameter  Error)  x 100  (1) 

A 

p 

where  p is  the  parameter  estimate  and  a is  the  estimate  of  the  la  confidence 
bound.  Thus,  EPE  is  the  estimated  parameter  error  in  percent  and  is  based  on  a 2a 
confidence  bound.  Experimental  results  at  NAVAIRTESTCEN  indicate  that  the 
range  of  EPE  values  presented  in  table  I can  be  used  to  express  qualitatively  the 
relative  accuracy  of  parameter  estimates. 

Table  I 


Qualitative  Assessment  of  Parameter  Estimate  Accuracy 


EPE 

(%) 

Qualitative  Parameter 
Estimate  Accuracy 

0 - 9.0 

Excellent 

9.0  - 20.0 

Good 

20.0  - 36.6 

Satisfactory 

36.6  and  above 

Unacceptable 

6)  NAVAIRTESTCEN  Technical  Report  SA-27R-76,  Verification  of  S-3A  Lateral- 
Directional  Power  Approach  Characteristics  Using  a Maximum  Likelihood 
Parameter  Identification  Technique,  Final  Report,  of  28  May  1976 
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RESULTS  AND  DISCUSSION 

DETERMINATION  OF  SHORT  PERIOD  RESPONSE  CHARACTERISTICS 


9.  The  aerodynamic  coefficients  describing  the  short  period  mode  of  motion  at 
constcmt  airspeed  (2a»  Mot»  Mq,  emd  were  determined  from  analysis  of 

pitch  doublet  cmd  sinusoidal  stick  input  responses.  Initially,  attempts  were  made  at 
including  the  veu'iation  of  normal  force  with  pitch  rate  (Zq)  in  the  estimation 
procedure.  However,  due  to  the  smeill  effect  of  this  parameter  on  the  total 
response,  the  estimated  Zq  exhibited  a very  low  confidence  factor.  Therefore,  Zq 
was  held  fixed  at  the  wind  tunnel  value  in  all  subsequent  estimations.  The  other 
aerodynamic  parameters  (those  that  effect  the  phugoid  mode)  in  the  equations  of 
motion  were  held  at  zero,  since  they  did  not  affect  the  short  period  response. 


10.  The  angle  of  attack  response  obtained  from  the  initial  estimations  exhibited  a 
Icig  of  approximately  0.2  sec  when  compared  to  the  measured  angle  of  attack 
response.  In  order  to  account  for  this  difference,  the  aircraft  equations  of  motion 
were  augmented  with  a description  of  the  measurement  lags  in  the  noseboom 
system.  This  resulted  in  a sensed  angle  of  attack  described  by  the  state  equation: 


1 1 \ 

. 1 

1 To.  ) 

“s  1 

^ Ta  } 

«t 

1 Ta  j 1 u ) 

with  the  measured  angle  of  attack  being  expressed  as: 


q (2) 


•m 


K. 


(3) 


In  this  manner,  the  lag  time  constant  in  the  measured  angle  of  attack  (1/Tqi)  was 
included  as  a parameter  in  the  estimation  procedure  and  resulted  in  better 
correlation  between  the  estimated  and  measured  angle  of  attack  responses. 

11.  The  nondimensional  aerodynamic  coefficients  obtained  from  the  maximum 
likelihood  analysis  of  the  short  period  responses  are  presented  in  appendix  D, 
figures  1 through  3.  The  range  of  expected  pcu-ameter  errors  resulting  from  the 
analysis  of  the  individual  maneuvers  is  presented  in  table  H.  The  identification 
procedure  yielded  excellent  accuracy  in  all  of  the  parameters  identified  with  the 
exception  of  Zg^  and  l/Ta  which  were  qualitatively  cissessed  as  good.  A comparsion 
with  the  wind  tunnel  results  of  reference  7 shows  excellent  agreement  in  the 
normed  force  coefficients  but  large  discrepcmcies  in  the  pitching  moment 
characteristics.  The  identified  results  indicate  that  the  airplane  is  less  statically 
stable  and  more  dynamically  stable  than  predicted  by  the  wind  tunnel  data.  In 
addition,  the  estimated  Mj  ^ shows  that  elevator  deflections  provide  a significantly 
larger  pitching  moment  than  that  estimated  by  the  wind  tunnel.  Typical  time 
history  compsu-isons  of  the  measured  aircraft  response,  the  maximum  likelihood 
estimated  response,  and  the  wind  tunnel  predicted  response  are  presented  in 
appendix  D,  figures  4 and  5.  The  maximum  likelihood  estimated  response  compares 
favorably  with  the  measured  flight  data  while  the  wind  tunnel  predicted  responses 
underestimate  the  magnitude  of  the  response  and  exhibit  less  damping. 


7)  GAC  Report  XA1128-104-1,  EA-6B  Basic  Stability  and  Control  Data  Report,  of 
February  1968. 
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Table  H 


Short  Period  Parameter  Estimate  Accuracy 


DETERMINATION  OF  LONG  PERIOD  RESPONSE 

12.  The  additional  derivatives  required  to  describe  the  long  period  mode  of  motion 
(Xjj,  Xj,,  Zy,  My)  were  determined  from  analysis  of  phugoid  responses.  Difficulty 
was  experienced  in  identifying  the  term  Xa  due  to  the  lack  of  cuigle  of  attack 
response  in  the  phugoid  mode.  However,  since  this  term  has  a significamt  impact  on 
the  determination  of  phugoid  damping  and  could  not  be  ignored,  it  was  calculated 
by  the  equation: 


and  held  fixed  in  the  identification  procedure.  The  other  terms  in  the  equations  of 
motion,  previously  determined  during  the  short  period  response  analysis,  were  held 
fixed  at  their  identified  value  during  the  long  period  identification  procedure. 


13.  The  phugoid  data  were  analyzed  over  a time  period  beginning  after  the  pilot's 
input  was  removed  (returned  to  the  trim  position)  and  continuing  through  one  and 
one-half  to  two  cycles  of  aircraft  motion.  Thus,  only  the  free  response  of  the 


phugoid  maneuver  was  used  in  the  analysis  (control  derivatives  X . , Z , and 


were  set  equal  to  zero).  The  initial  force  and  moment  conditions  (X^)  were 
estimated  at  the  begiiming  of  each  time  period  to  take  into  accound  the  fact  that 
the  analysis  ^as  initiated  from  an  off  trim  condition  (see  appendix  B for  a 
definition  of  X ).  The  nondimensional  aerodynamic  coefficients  obtained  from  the 
long  period  analysis  are  presented  in  table  HI.  The  identification  procedure  yielded 
excellent  results  in  the  estimation  accuracy  of  Cxy  and  and  good  results  for 
CMy.  ^ time  history  comparison  of  the  measured  flight  data  and  the  maximum 
likelihood  estimated  response,  presented  in  appendix  D,  figure  6,  shows  excellent 


agreement  in  all  of  the  measured  responses. 
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Table  HI 

Aerodynamic  Coefficients  Describing  Variation  in  Force  and 
Moment  Characteristics  with  Forward  Velocity 

Loading  - 3 TJ  Pods  and  2 Drop  Tanks 
True  Angle  of  Attack  - 14  deg 
Landing  Gear  - DOWN 
Flaps  - 30  deg 


Parameter 

Value 

EPE  (%) 

-.213 

2.9 

^Zu 

-1.615 

1.3 

CMu 

-.01 

18.7 

VERIFICATION  OF  IDENTIFIED  PARAMETERS 


14.  The  aerodynamic  coefficients  identified  from  the  maximum  likelihood  estima- 
tion were  verified  through  determination  of  the  accuracy  of  the  identified 
parameters  (as  based  on  estimated  parameter  error),  time  history  comparisons  of 
the  measured  and  estimated  responses,  and  compcu-ison  of  computed  response 
characteristics  with  conventional  flight  test  results.  The  parameter  accuracy  (EPE) 
and  time  history  comparisons  made  indicate  that  a high  level  of  confidence  should 
be  placed  on  the  identified  pcu-ameters.  A discussion  of  computed  versus  measured 
chcuracteristics  follows. 


15.  The  frequency  and  damping  characteristics  of  both  the  short  and  long  period 
responses  were  obtained  by  first  determining  the  Eigenvalues  of  the  identified  state 
matrix,  F,  and  then  calculating  undamped  natural  frequency  and  damping  ratio 
from  the  real  and  imaginary  parts  of  those  Eigenvalues.  The  resulting  short  period 
frequency  and  damping  characteristics  are  presented  in  appendix  D,  figures  7 cuid  8. 
A comparison  of  the  short  and  long  period  maximum  likelihood,  conventional  flight 
test  data  analysis,  and  wind  tunnel  estimated  frequency  and  damping  characteris- 
tics is  presented  in  table  IV.  Excellent  agreement  between  the  maximum  likelihood 
estimated  and  conventional  data  analysis  results  are  indicated  at  all  test 
conditions.  The  wind  tunnel  estimates  provide  a good  determination  of  the  short 
period  frequency  characteristics  but  considerably  underestimate  the  damping.  This 
result  is  to  be  expected  from  the  much  lower  pitch  damping  term  (Cj^^  + Cj^.) 

obtained  from  the  wind  tunnel  data  (see  paragraph  11).  *1  a 
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Table  IV 


Frequency  and  Damping  Characteristics 


Loading  - 3 TJ  Pods  and  2 Drop  Tamks 
Landing  Gear  - DOWN 
Flaps  - 30  deg 


True 
Angle  of 
Attack 
(deg) 

Modal  Parameter 

Maximum 

Likelihood 

Identification 

Result 

Conventional 

Data 

Analysis 

Result 

Wind 

Tunnel 

Estimate 

12 

Short 

Period 

Wrisp 

(rad/sec) 

1.05 

1.0 

1.05 

^SP 

0.66 

0.61 

0.45 

18.5 

Short 

Period 

U)nsp 

(rad/sec) 

0.88 

1.0 

1.04 

^SP 

0.67 

0.66 

0.34 

14 

Long 

Period 

"•np 

(rad/sec) 

0.148 

0.147 

- 

■■ 

0.092 

0.093 

- 

16.  Additional  verification  of  the  maximum  likelihood  estimated  parameters  was 
obtained  by  comparison  with  conventional  static  longitudinad  stability  test  results. 
Static  longitudinal  stability  characteristics  in  terms  of  elevator  position  gradients 
with  airspeed  can  be  determined  by  evaluating  the  forward  airspeed  to  horizontal 
stabilizer  transfer  function  at  steady  state  conditions: 


u 

_ 

K {=  * ''"i)  (*  * 

(5) 

A 

(s^  + 2Cp  (Drip  s + '^np^)(^^  + 2^SP  '^^SP 

S + 'Dpgp  j 

Evaluating  this  equation  at  s = o (steady  state)  yields 


u(s) 

^e(s) 


“isp^ 


(6) 
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which  is  equivalent  to 


in  terms  of  the  identified  pcirameters.  Equation  (7)  can  be  interpreted  as  the 
gradient  of  airspeed  with  elevator  position  during  a static  longitudinal  stability 
test.  The  static  longitudinal  stability  results  utilizing  the  estimated  peu-ameters 
obtained  from  dynamic  stability  tests  are  compared  with  the  conventionad  results 
of  static  longitudinal  tests  (reference  8)  in  appendix  D,  figure  9.  Again,  excellent 
agreement  is  obtained  between  the  two  methods,  thereby  lending  further  credence 
to  the  parameters  identified  by  the  maximum  likelihood  algorithm. 


8)  NAVAIRTESTCEN  Technical  Report  FT-68R-73,  Flying  Qualities  and  Perform- 
ance Technical  Evaluation  of  the  EA-6B  Airplame,  Final  Report,  of  30  October 
1973 


9 


TM  78-2  SA 


CONCLUDING  REMARKS 

17.  The  stability  and  control  derivatives  describing  the  longitudinal  aerodynamics 
of  the  EA-6B  airplane  in  the  catapult  launch  configuration  have  been  accurately 
determined  from  flight  data  through  the  use  of  a maximum  likelihood  estimation 
algorithm.  The  identified  derivatives  were  obtained  from  a trimmed  level  flight 
condition  and  are  useful  for  point  simulations  and  specification  determination.  If 
the  simulation  of  out  of  trim  flight  conditions  is  required,  tests  would  have  to  be 
conducted  at  other  them  level  flight  conditions  in  order  to  determine  the  effects  of 
thrust  and  horizontal  stabilizer  trim  settings  on  the  aerodynamic  parameters. 

18.  The  aerodynamic  coefficients  were  generated  in  a form  suitable  for  the  i 

construction  of  a data  base  for  simulation  of  EA-6B  airplane  characteristics.  This  { 

data  biise  has  direct  application  to  a parametric  study  of  catapult  laimch  i 

capabilities  of  the  EA-6B  airplane  through  computer  simulations.  The  data  base  can 

also  be  used  for  optimization  of  control  system  configurations  in  ACLS  simulations,  i 

for  verifying  or  constructing  the  data  base  of  the  EA-6B  Operational  Flight 

Trainer,  and  for  determining  aircraft  modal  peirameters  for  specification  analysis.  | 

j 
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MOMENT  OF  INERTA  ABOUT  THE  PITCH  AXIS 
THOUSAND  SLUG -FT^ 
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EQUATIONS  OF  MOTION 
LINEAR  STATE  SPACE  MODEL 

1.  The  longitudinal  equations  of  motion  and  measurement  equations  used  in  the 
maximum  likelihood  parameter  identification  algorithm  were  written  in  steuidaurd 
linear  state  space  format  in  a body  fixed  axis  system. 

Equation  of  Motion; 

X = Fx  + Gu  + X 

o 

where: 

X is  an  nxl  state  vector 
u is  a pxl  control  vector 
F is  a nxn  matrix  of  stability  derivatives 
G is  a nxp  matrix  of  control  derivatives 

X is  an  nxl  vector  of  initial  conditions  on  the  state  vector 
o 

Measurement  Equation: 

y = Hx  + Du  + B + V 
where: 

y is  an  rxl  measurement  vector 
b is  an  rxl  measurement  bias  vector 
V is  an  rxl  measurement  error  vector 

H is  an  rxn  matrix  relating  the  measurement  to  the  state  vector 
and: 

n is  the  number  of  states 
p is  the  number  of  controls 
r is  the  number  of  measurements 

This  type  of  formulation  is  very  general  and  2dlows  for  the  maximum  amount  of 
flexibility  in  including  the  correct  relationship  between  the  dynamic  equations  of 
motion  and  measurements  of  aircraft  response  v2U‘iables.  Thus,  errors  in  instrumen- 
tation system  measurements  are  included  as  an  integral  part  of  the  identification 
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process  (measurement  errors  are  estimated  simultaneously  with  aircraft  stability 
and  control  derivatives).  For  example,  during  the  identification  process,  angle  of 
attack  measurements  are  corrected  for  pitch  rate  effects  on  the  local  flow  field  at 
the  nose  boom  emd  instrument  bicis,  measurement  random  error,  and  angle  of  attack 
vane  scale  factor  are  estimated. 


LONGITUDINAL  MODEL 

2.  For  the  longitudinal  case,  the  state  space  model  takes  on  the  following  form: 


Equations  of  Motion: 


State  vector:  x = |^a  , u,  q,  0 , j 


where  a , the  sensed  angle  of  attack  was  included  to  account  for  lags  in  the  angle 
of  attaci  sensing  system  of  the  test  airplane. 


Control  vector:  u “ 
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TEST  INSTRUMENTATION  PARAMETERS 
EA-6B  Airplane 
BuNo  156478 


Parameter  Cockpit 

Nose  Boom  Airspeed 
Nose  Boom  Altitude 
Nose  Boom  Angle  of  Attack 


Nose  Boom  Sideslip  Angle 

Production  Airspeed  X 

Production  Altitude  X 

Production  Angle  of  Attack  X 

Outside  Air  Temperature  X 


Left  and  Right  Turbine  Discharge  Pressure 
Left  and  Right  Low  Speed  Compressor  RPM 
Left  and  Right  High  Speed  Compressor  RPM 
Left  and  Right  Turbine  Outlet  Temperature 
Left  and  Right  Fuel  Used  Counters  (Volume) 
Horizontal  Stabilizer  Position 
Pitch  Angle 

Longitudinal  Stick  Position 
Pitch  Rate 

CG  Normal  Acceleration 
Flaperon  Position 
Time  Code 


Mag  Tape 
and  TM 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 
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(BAD  ')  (bad  ’)  (BAD  ')  (BAD  ')  (BAD ’) 
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LOADING  - 3 TJ  POOS  AND  2 OBOP  TANKS 
LANDING  GEAB  • DOWN 
FLAPS -30  DEG 
SAS-OFF 

ALTITUDE  - 7,000  FT  (2  130  m) 

C6  POSITION -24.4XIHAC 
O OOUBL'ET  INPUTS 
□ SINE4IIAVE  INPUTS 


POWEB  OFF  LIFT  COEFFICIENT 


Figure  1 
EA-6B  Airplane 
BuNo  156478 

Short  Period  Aerodynamic  Coefficients 
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LANDING  GEAR  - DOWN 
FLAPS -30  DEG 

□ O 3TJ  PODS  AND  2 DROP  TANKS 
• CLEAN  LOADING 


POWER  OFF  LIFT  COEFFICIENT 


Figure  3 
EA-6B  Airplane 
BuNo  156478 


Angle  of  Attack  and  Stabilizer  Characteristics 


ELEVATOR  NORMAL  MTCH  PITCH  AH6LE  OF 

POSITION  ACCELERATION  ATTITUDE  HATE  ATTACK 

RAO  FT/SEC**2  RAO  RAO/SEC  RAO 
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LOADING  - 3 TJ  PODS  AND  2 DROP  TANKS 
LANDING  GEAR  - DOWN 
FLAPS  -30  DEG 
SAS • OFF 

ALTITUDE  - 7,000  FT  (2  130  m) 

CG  RANGE  - 23.6  TO  25.2%  MAC 

GROSS  WEIGHT  RANGE  - 46,500  TO  45,340  LB  (21  090  TO  20  570  kg) 


□ O NATC  MAXIMUM  LIKELIHOOD  ESTIMATE 
• CONVENTIONAL  DATA  ANALYSIS 
A WIND  TUNNEL  PREDICTION 


POWER  OFF  LIFT  COEFFICIENT 


‘ Figure  7 

EA-6B  Airplane 
BuNo  156478 

Short  Period  Frequency  and  Damping  Characteristics 
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DAMPING  RATIO  UNDAMPED  NATURAL  FREQUENCY 


nai-oaa 

NOlilSOd 

uazmavis 


Figure  9 
EA-6B  Airplane 
BuNo  156478 

Comparison  of  Static  Longitudinal  Stability  Characteristics 
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UST  OF  SYMBOLS 


Symbol 

Definition 

Units 

a 

z 

Vertical  acceleration 

ft/sec^ 

(m/s^ 

b 

Measurement  bias 

c 

Mean  aerodynamic  chord 

ft 

(m) 

X 6 

Nondimensional  partied  derivative  of 

K force  or  moment  (X,  Z,  M)  with 
respect  to  state  vcuriable,  x,  or  control, 

6 

PWR  OFF 

Power-off  lift  coefficient 

- 

- 

D 

Matrix  relating  measurements  to 
control  vector 

EPE 

Estimated  Parameter  Error 

% 

F 

Matrix  of  stability  derivatives 

G 

Matrix  of  control  derivatives 

g 

Acceleration  due  to  gravity 

ft/sec^ 

(m/s^! 

H 

Matrix  relating  measurements  to 
state  vector 

lyy 

Moment  of  inertia  about  pitch  axis 

slug-ft^ 

(kg-m 

Ka 

Angle  of  attack  vane  scale  factor 

1 

X 

X distance  of  angle  of  attack  vane 
to  eg 

ft 

(m) 

Iz 

Z distance  of  angle  of  attack  vcuie 
to  eg 

ft 

(m) 

M 

Pitching  moment  about  Y axis 

ft-lb 

(N-m) 

Partial  derivative  of  pitching  moment 
with  respect  to  state  variable,  x,  or 
control,  6 

A 

P 

Parameter  estimate  (general) 

q 

Pitch  rate 

rad/sec 
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Symbol 

Definition 

Units 

q 

Dynamic  pressure 

Ib/ft^ 

(Pa) 

s 

Laplace  operator 

- 

S 

Wing  area 

ft^ 

(m^) 

t 

Time 

sec 

u 

Longitudinal  component  of  velocity 

ft/sec 

(m/s) 

V 

Free  stream  velocity 

ft/sec 

(m/s) 

w 

Vertical  component  of  velocity 

ft/sec 

(m/s) 

X 

State  variable  ( a , u,  q,  0, 

- 

X 

X component  of  force 

lb 

(N) 

X , 

x’ 

^6 

Partial  derivative  of  X component  of 
force  with  respect  to  state  variable, 

X,  or  control,  6 

y 

Measurement  variable  (a  , u , a , 

_ , “ ™ “ 

®m’  ' 

m 

z 

Z component  of  force 

lb 

(N) 

^6 

Partied  derivative  of  Z component  of 
force  with  respect  to  state  variable, 

X , or  control , 6 

(') 

Time  rate  of  ch£mge 

sec~^ 

A 

( ) 

Estimate 

(■) 

Vector 

rad 

a 

Angle  of  attack 

6 

Control  deflection 

rad 

i 

s 

Horizontal  stabilizer  position 

rad , deg 

C 

Damping  ratio 

- 

0 

Pitch  attitude 

rad 

o 

Parameter  estimate  v2U'iance  or  con- 
fidence bound 

- 

T 


Time  constant 
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Symbol  Definition  Units 

V Measurement  random  error 

Natrual  frequency  rad/sec 


Measured 

True 

Trim  condition 

Inverse 

Transpose 
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